Dynamical mean-field approximations for a diffusive pair contact process 
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Dynamical mean-field approximations are performed to study the phase transition of a pair contact process 
with diffusion in different spatial dimensions. The level of approximation is extended up to 18-site clusters for 
the one-dimensional model. The application of coherent anomaly method shows that the critical /3 exponent 
does not depend on the strength of diffusion rate. The extension of the mean-field approximation to higher 
dimensions also suggests that the critical behavior may be described by a unique set of exponents. 
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The study of nonequilibrium critical phenomena has re- 
ceived considerable attention but their understanding is far 
from complete. Even the simplest models can be theoreti- 
cally challenging. A well-known example is the pair contact 
process model with diffusion (PCPD) |1]. Despite its inten- 
sive study there is no satisfying consensus about the accu- 
rate values of critical exponents (see Ref. |2|, for a recent 
review and further references). Different methods and ap- 
proximations were applied, such as Monte Carlo simulation 
(MC) 1 3, 4,'?','^, dynamical mean-field approximation (DMF) 
11, 6], density-matrix renormalization group study | IJ, and 
field-theoretic approach Most of the studies have been 

devoted to the one-dimensional case but the published numer- 
ical results are partially contradictory. In addition to the sim- 
ulation difficulties, a consistent field theory description of the 
model is still an unresolved issue |8]. Instead of summing 
all the raised proposals |2], we consider one of the sugges- 
tions that the value of critical exponent may depend on the 
strength of diffusion. This idea was inspired by different MC 
simulations and DMF approximations 10]. This latter 

approach, however, left some ambiguities for strong or very 
weak diffusion rate [6]. Subsequent simulations I IQifTlll re- 
vealed the concept of a novel universality class with a unique 
set of exponents. However exact calculation, albeit for a 
slightly different model, suggested the possibility of two tran- 
sitions depending on the diffusion rate in higher spatial dimen- 
sions 1 12]. The primary purpose of this Brief Report to clarify 
the mentioned ambiguities of DMF approach by significantly 
improving the level of one-dimensional approximations and 
to extend the calculations to higher spatial dimensions. 

The simplicity of PCPD model, especially in one dimen- 
sion, makes possible to extend the DMF approximation to 
large-size clusters. Such a large number of approximations 
for the order parameter function enables the application of the 
coherent anomaly method (CAM) to derive quantitative pre- 
diction to the value of [3 exponent ] 13]. Another challenge is 
to develop the DMF approximation to higher spatial dimen- 
sions. Generally, such kind of calculation is expected to yield 
qualitatively correct results in high dimensions, therefore it 
can help us to clarify the raised question whether the values 
of critical exponents depend on the strength of diffusion rate. 



The PCPD model J 1 , 14j is the extension of the pair-contact 
process model (PCP) \15\ that belongs to the well-known DP 
universality class lllal . The PCPD model consists of a d- 
dimensional lattice, with periodic boundaries, where each site 



FIG. 1: Dynamical mean-field approximations of the particle con- 
centration at D = 0.1 forn = 2, 4, 6, 18-point levels (from right 
to left curves). The arrow shows the location of the extrapolated 
phase transition point. 



is either vacant or occupied by a single particle. The dynam- 
ical rules for updating the system are defined as follows. A 
randomly chosen pair of nearest-neighbor particles is annihi- 
lated with a probability p(l — D). An additional particle is 
created around the given pair with probability (1 — — D) 
if it is not forbidden by the exclusion principle. Furthermore 
individual particles are allowed to hop to the nearest-neighbor 
empty site with diffusion rate D. This version is also called 
restricted PCPD model because the exclusion law prohibits 
double occupancy. In this model the order parameter is the 
concentration of pairs of particles (u) that becomes zero above 
a critical value of p when the system arrives to one of the two 
possible absorbing states 0]. 

The application of DMF approximation, which is a dynam- 
ical version of the cluster variation method, involves finding 
a hierarchy of evolution equations for the probability distri- 
butions of configurations within a cluster of n sites (for de- 
tails, see e.g. Refs. ] 17, 18]). The calculation of transitions of 
cluster probabilities would generate infinite hierarchy of equa- 
tions, since transitions in an n-site cluster may depend on sites 
outside the cluster To avoid this the n-site approximation es- 
timates m-site probabilities (for to > n) by using Bayesian 
extension process. This type of approximation has proved to 
be successful in the evaluation of phase diagram for many dif- 
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TABLE I: Results of n-site approximations for the one-dimensional model. The extrapolated values of Pc and the order parameter 
exponents p for different values of diffusion rate (D) are also given. Numbers in parentheses correspond to the uncertainty in the 
last digit. 
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4 


0.157167(1) 


48 4(1) 


0.160964(1) 


82.1(2) 


0.209694(2) 


51.7(1) 


0.294937(1) 


17.61(1) 


6 


0.137271(1) 


94.2(4) 


0.140885(1) 


147.3(6) 


0.184175(4) 


93.2(2) 


0.274965(2) 


17.40(2) 


8 


0.129264(2) 


127.1(6) 


0.132484(4) 


207.6(10) 


0.171820(6) 


139.2(8) 


0.260674(6) 


19.63(2) 


10 


0.124488(5) 


158.2(10) 


0.127542(4) 


265.5(18) 


0.164401(6) 


190.9(10) 


0.250095(8) 


23.08(3) 


12 


0.121299(7) 


205.9(20) 


0.124240(10) 330.8(16) 


0.159422(8) 


248.9(14) 


0.242122(12) 


26.83(4) 


14 


0.119015(7) 


254.8(20) 


0.121799(10) 413.6(30) 


0.155843(10) 308.0(20) 


0.235710(16) 


33.31(10) 


16 


0.11730(3) 


308.2(60) 


0.11996(3) 


489.4(60) 


0.15309(4) 


367.2(40) 


0.23071(3) 


41.02(20) 


18 


0.11596(4) 


370.5(50) 


0.11841(11) 


590.2(80) 


0.15053(10) 


453.3(70) 


0.22676(12) 


49.19(40) 


Pc 


0.1053(2) 




0.1066(5) 




0.1335(1) 




0.1955(5) 




p 




0.54(4) 




0.55(6) 




0.51(4) 




0.48(3) 



ferent nonequilibrium models (TlSIIliSlilllllllll^l. 

Evidently, this method was also applied to the present model 
in one dimension up to n = 7, however, the corresponding 
results were insufficient to derive reliable prediction to the 
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values of critical exponents ^ [Till . Furthermore CAM re- 
sult could not be presented for strong {D > 0.5) and for weak 
(D < 0.1) diffusion rates fl. 



To eliminate these shortages now the approximation is ex- 
tended up to n = 18 level. Naturally, for higher level approx- 
imations the equations can be only build and solved numeri- 
cally 1 25, 27]. The stationary distribution is attained via nu- 
merical integration when the sum of the absolute values of all 
time derivatives smaller than e (generally e — 10^^ is used). 
Increasing the level of approximation the necessary number 
of iterations to attain the steady-state configurations probabil- 
ities increases drastically. To demonstrate the speed of conver- 
gence for large n, four-week running is necessary to get data 
at 18-point level. Despite the time consuming convergence 
during the iterations, the calculations are tractable even by a 
personal computer. Figure [l] illustrates the improvement of 
solutions as we increase the levels of approximation. Here the 
particle concentration functions p{p) are plotted for D — 0.1, 
which show linear decay around the transition point. The con- 
centration of pairs can be displayed in a similar figure but u{p) 
decreases quadratic ally around p" for n > 2. 

The order-parameter exponent can be estimated from the 
mean-field data following the CAM approach developed by 
Suzuki [13.1 . In the vicinity of the critical point the order pa- 
rameter function can be fitted as u{p) « a„ (p" — p)'^ , 
where p" denotes the value of the critical point at 7i -point level 
and Pc is the extrapolated value to the n ^ oo limit. As we 
mentioned /3*^^ = 2 is valid for all values of D if n > 2. Ta- 
ble|l]summarizes the results for the estimated transition points 
(p") and amplitudes (a„) at different levels of approximation 
for various values of diffusion rate. (The starting point of fit- 
ting is determined when the corresponding particle concentra- 
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FIG. 2: CAM scaling of the critical mean-field coefficients for the 
order parameter at different values of the diffusion constant. The 
lines correspond to the — 0.5 exponent. The results of n = 
4, 6, 8, 18-point level of approximations are plotted from right to 
left. 

tion drops below 0.01.) According to the CAM analysis the 
amplitudes a„ scales as 

^ (p: - Pcf-^-r (1) 

As a consequence, the log — log plot of the amplitudes vs. 
Pc ^ Pc yields an estimation for the (3 exponent. Figure |2] 
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shows these plots for different values of D. This figure indi- 
cates clearly that wrong value of /3 could be deduced by this 
technique if the calculations are limited to small-size clusters. 
For large n, however, the power law fits of the correspond- 
ing amplitudes become possible. The deviation of a„ from 
an expected power law function is pronounced for large D 
when n are small. Although the convergence of a„ is ob- 
vious, still the uncertainty of the estimated (3 value is close 
to 10%. The extrapolated values of pc and the [3 exponents 
for different diffusion rates are also given in Table |I] Within 
the mentioned uncertainty the same f3 = 0.50(5) exponents 
can be fitted for all values of D as displayed in Fig.|2] This 
P value is comparable to some earlier published MC results 
0,1^. To end this section it is worth emphasizing the sensitiv- 
ity of CAM on the levels of approximation. The comparison 
of Fig. n and Fig. |2] underlines the different behaviors of 
and a„ as a function of n. As an example, the approach of a„ 



to the power law function is significant if we increase n from 
8 to 18 for D — Q.l but the location of the transition point 
changes less then 11%. Similar behavior was reported by a 
recent study f28] where the same model was investigated by 
DMF approximation up to n = 13 for D — 0.5. 

In the light of a field-theoretic calculation fT2il it is also in- 
teresting to extend our investigations to higher spatial dimen- 
sions (d > 1). The main question is whether two transitions 
exist as a function of diffusion rate, D. The lowest level that 
considers the value of d is the pair approximation {n = 2). In 
spite of its failure in one dimension 1 1, 2], this approximation 
is expected to provide qualitatively correct behavior in higher 
dimensions. At this level the two independent variables are 
the particle concentration denoted by p and the concentration 
of pairs u. The equations of their time evolution on a d di- 
mensional lattice can be written as 



-2{1 - D)pu + {1 - D){1 - p) 



u{p — u) 



(2) 



+ {l-D){l-p) 



u(p — u) 



k=0 



{u)Hp~u) 



2m — A; 



Y.{k+i)r){p-uf{i-2p 



\m — k 



ra / \ / m 

+-;4^i:(";)(.-")'u-v+.r-'(i:(/-.) 



u\p-uy 



(3) 



r 



where m = 2(i — 1. By comparing these equations in case of 
d = 1 to the equations in Ref. |2'l, one can observe that a pref- 
actor 2 is missing from the diffusion part of Eq.|3] The slight 
difference can be explained as follows. When doing MC sim- 
ulation an elementary diffusion process includes to choose a 
particle and a direction as well. Afterwards the jump is ex- 
ecuted with probability D if the target site is empty. There- 
fore, the probability to jump to a specified direction is equal 
to D/ {2d). The normalization is fulfilled by choosing the lat- 
ter jumping rate. 

Returning now to the above equations the solutions in the 
active phase are 



1-p 
l-Sp 



where 



(4) 



l~Sp\ 2m- {I - D){2m3p+l -p) 
1-p J 2m - (1 - D){2m 3p + 1 - 3p) 



As a consequence, the value of critical point is 



pt'^{d,D) = 



4rf-3+D 
(12d-7)(l-D) 



if 0< D < 



otherwise . 



2(3d-l) 



(5) 



It means that, as in one dimension, two different regions can 
be detected at any finite values of d as a function of D. Nev- 
ertheless, the region that characterizes the high-Z? behavior, 
including pc = 1/3, expands if we increase d and in the infi- 
nite dimension limit the solution convergences to the solution 
of the one-site approximation. From this result one may con- 
jecture that there is no different universality classes at low and 
high diffusion rate for the restricted PCPD in higher dimen- 
sions. 

More convincing arguments can be achieved by performing 
DMF approximation for larger clusters in two dimensions. In 
the absence of diffusion (D — 0) the two-point approximation 
yields qualitatively incorrect result, since the particle concen- 
tration becomes zero at the transition point (p^P — 5/17 w 
0.2941). This is contrary to the MC result that revealed a fi- 
nite 'natural' density, pnat = 0.1477 at the transition point 
{pc = 0.2005) for the PCP model f^^l. The four-site approx- 
imation (on 2x2 clusters) yields a better estimation for the 
transition point (p'^P = 0.2691(3)), however, it also predicts 
zero particle concentration at the transition point. This latter 
failure is already diminished at nine-site (3x3 clusters) ap- 
proximation. As the inset of Fig.|3]shows, this level of approx- 
imation can already describe qualitatively well the behavior of 
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FIG. 3: The pair-density function of the two-dimensional model at 
9-point level for D = 0, 0.1, 0.5 and 0.9 values (from left to right 
curves). The inset shows particle density at D = forn = 2, 4, and 
9-point levels (from right to left curves). 

the two-dimensional PCP model: it gives pnat — 0.1197(2) 
at p^P — 0.22607(3). The reason why the 9-point level gives 
a qualitatively different result may be explained by the fact 
that the 3x3 cluster size is the smallest cluster that covers 
all the elementary processes. (For example, a particle creation 
at the ends of an occupied pair "requires" a minimum size of 
3 in one direction.) Similar qualitative improvement of the 
results have already been observed in a two-offspring branch- 
ing annihilating random walk model |30]. In the latter one- 
dimensional model at least 5-point level of approximation was 
necessary for the correct description. The success of 9-point 



level of approximation also means that the smallest cluster 
that may give correct description of our PCPD model in three 
dimensions consists 3 x 3 x 3 = 27 lattice points. This level 
is not obtainable presently. Returning to the two-dimensional 
case, the 9-point approximation does not predict qualitatively 
different behavior for low and high diffusion rate. Figure |3l 
demonstrates that the value of critical point increases with D 
{plP{D = 0.1) 0.2607(5), pIp{D 0.9) = 0.3315(5)) 
while the order parameter function does not differ qualita- 
tively for _D = 0.1 and D — 0.9. This result also supports 
the single universality class conjecture. 

To summarize, DMF approximation is extended to large- 
cluster levels to study the critical behavior of the one- 
dimensional restricted PCPD model. We find that the ambigu- 
ity of CAM approach for weak and strong D can be resolved 
by increasing the level of approximation. The approximation 
suggests uniform f3 exponent for all values of D within the un- 
certainty of the estimations. This /? = 0.50(5) value is com- 
parable with some published MC results and supports that the 
model leaves the DP universality class and described by an 
exponent that also differs from the PC class. The approxima- 
tions in higher dimensions also support the conjecture that the 
PCPD model can be described by a single universality class 
that is independently from the rate of diffusion. 

Acknowledgments 

I thank Geza (Ddor (who suggested this investigation), and 
Gyorgy Szabo for illuminating discussions. This work was 
supported by the Hungarian National Research Fund under 
Grant No. T-47003. 



[1] E. Carlon, M. Henkel, and U. Schollwock, Phys. Rev. E 63, 
036101 (2001). 

[2] M. Henkel and H. Hinrichsen, J. Phys. A 37, Rl 17 (2004). 

[3] H. Hinrichsen, Phys. Rev. E 63, 036102 (2001). 

[4] M. Henkel and U. Schollwock, J. Phys. A 34, 3333 (2001). 

[5] K. Park and I.-M. Kim, Phys. Rev. E 66, 027106 (2002). 

[6] G. Odor, Phys. Rev. E 62, R3027 (2000). 

[7] F. van Wijland, U. C. Tauber, and O. Deloubriere, cond- 

mat/03 11568 (2003). 
[8] H.-K. Janssen, F. van Wijland, O. Deloubriere, and U. C. 

Tauber, condmat/0408064 (2004). 
[9] R. Dickman and M. A. F. de Menezes, Phys. Rev. E 66, 045101 

(2002) . 

[10] J. Kockelkoren and H. Chate, Phys. Rev. Lett. 90, 125701 

(2003) . 

[11] G. Odor, Phys. Rev. E 67, 0161 1 1 (2003). 
[12] M. Paessens and G. M. Schutz, J. Phys. A 37, 4709 (2004). 
[13] M. Suzuki, J. Phys. Soc. Jpn. 55, 4205 (1986). 
[14] M. J. Howard and U. C. Tauber, J. Phys. A 30, 7721 (1997). 
[15] I. Jensen, Phys. Rev. Lett. 70, 1465 (1993). 
[16] M. A. Mufioz, G. Grinstein, R. Dickman, and R. Livi, Phys. 
Rev. Lett. 76, 451 (1996). 



[17] R. Dickman, Phys. Lett. A 122, 463 (1987). 

[18] J. Marro and R. Dickman, Nonequilibrium Phase Transitions 

in Lattice Models (Cambridge University Press, Cambridge, 

1999). 

[19] H. A. Gutowitz, J. D. Victor, and B. W. Knight, Physica D 28, 
18 (1987). 

[20] G. Szabo, A. Szolnoki, and L. Bodocs, Phys. Rev. A 44, 6375 
(1990). 

[21] D. ben Avraham and J. Kohler, Phys. Rev. A 45, 8358 (1992). 

[22] G. Szabo and I. Borsos, Phys. Rev. E 49, 5900 (1994). 

[23] M. C. Marques, M. A. Santos, and J. F. F. Mendes, Phys. Rev. 

E65, 016111 (2002). 
[24] A. Szolnoki, Phys. Rev. E 62, 7466 (2000). 
[25] A. Szolnoki, Phys. Rev. E 66, 057102 (2002). 
[26] G. Szabo and C. Toke, Phys. Rev. E 58, 69 (1998). 
[27] R. Dickman, Phys. Rev. E 66, 036122 (2002). 
[28] S.-C. Park and H. Park, condmat/04091 15 (2004). 
[29] J. K. L. da Silva and R. Dickman, Phys. Rev. E 60, 5 126 (1999). 
[30] D. Zhong, D. ben Avraham, and M. A. Muiioz, Europ. Phys. J. 

B 35, 505 (2003). 



